Contribution of host species and pathogen clade to snake fungal disease hotspots in Europe.

Infectious diseases are influenced by interactions between host and pathogen, and the number of infected hosts is rarely homogenous across the landscape. Areas with elevated pathogen prevalence can maintain a high force of infection and may indicate areas with disease impacts on host populations. However, isolating the ecological processes that result in increases in infection prevalence and intensity remains a challenge. Here we elucidate the contribution of pathogen clade and host species in disease hotspots caused by Ophidiomyces ophidiicola, the pathogen responsible for snake fungal disease, in 21 species of snakes infected with multiple pathogen strains across 10 countries in Europe. We found isolated areas of disease hotspots in a landscape where infections were otherwise low. O. ophidiicola clade had important effects on transmission, and areas with multiple pathogen clades had higher host infection prevalence. Snake species further influenced infection, with most positive detections coming from species within the Natrix genus. Our results suggest that both host and pathogen identity are essential components contributing to increased pathogen prevalence.

Isolation of the entomopathogenic fungus Beauveria bassiana from a skin lesion in a wild Nikolsky's viper (Vipera berus nikolskii).

Mycotic lesions of the skin of snakes are often associated with the fungus Ophidiomyces ophiodiicola, but other pathogens can cause similar signs. A skin sample from a wild Nikolsky's viper (Vipera berus nikolskii) with dermal lesions was collected in eastern Ukraine. A pure fungal culture was obtained and identified using nucleotide sequence analysis as the entomopathogenic species Beauveria bassiana sensu lato. Although Beauveria spp are considered to be non-pathogenic in vertebrates, sporadic infections have been reported. This report highlights the need to explore various pathogens when diagnosing the cause of snake integumentary lesions.

Cold Finger: Raynaud Phenomenon Following Snakebite Envenoming by Nikolsky's Viper (Vipera berus nikolskii).

A field biologist was bitten by a female Nikolsky's viper (Vipera berus nikolskii) in Kharkiv Oblast, Ukraine. Two months later, the patient began to experience cold-induced vasospasm of the affected digit diagnosed as acquired Raynaud phenomenon. The patient had more than 30 occurrences during the single winter following the bite, but the signs and symptoms of Raynaud phenomenon disappeared with the end of winter. This report describes the case and puts it into context with the literature on the topic of toxin-induced peripheral vasospastic disorders and their potential importance in snakebite envenoming.

Evolutionary allometry and ecological correlates of fang length evolution in vipers.

Traits for prey acquisition form the phenotypic interface of predator-prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.

Genetic Diversity of Porcine Circovirus 2 in Wild Boar and Domestic Pigs in Ukraine.

Porcine circovirus type 2 (PCV2) is responsible for a number of porcine circovirus-associated diseases (PCVAD) that can severely impact domestic pig herds. For a non-enveloped virus with a small genome (1.7 kb ssDNA), PCV2 is remarkably diverse, with eight genotypes (a-h). New genotypes of PCV2 can spread through the migration of wild boar, which are thought to infect domestic pigs and spread further through the domestic pig trade. Despite a large swine population, the diversity of PCV2 genotypes in Ukraine has been under-sampled, with few PCV2 genome sequences reported in the past decade. To gain a deeper understanding of PCV2 genotype diversity in Ukraine, samples of blood serum were collected from wild boars (n = 107) that were hunted in Ukraine during the November-December 2012 hunting season. We found 34/107 (31.8%) prevalence of PCV2 by diagnostic PCR. For domestic pigs, liver samples (n = 16) were collected from a commercial market near Kharkiv in 2019, of which 6 out of 16 (37%) samples were positive for PCV2. We sequenced the genotyping locus ORF2, a gene encoding the PCV2 viral capsid (Cap), for 11 wild boar and six domestic pig samples in Ukraine using an Oxford Nanopore MinION device. Of 17 samples with resolved genotypes, the PCV2 genotype b was the most common in wild boar samples (10 out of 11, 91%), while the domestic pigs were infected with genotypes b and d. We also detected genotype b/d and b/a co-infections in wild boars and domestic pigs, respectively, and for the first time in Ukraine we detected genotype f in a wild boar from Poltava. Building a maximum-likelihood phylogeny, we identified a sublineage of PCV2 genotype b infections in both wild and domestic swine, suggesting a possible epizootic cluster and an ecological interaction between wild boar and domestic pig populations in northeastern Ukraine.

Evolutionary allometry and ecological correlates of fang length evolution in vipers

Traits for prey acquisition form the phenotypic interface of predator–prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems.

Características para aquisição de presas formam a interface fenotípica das interações predador-presa. Em predadores venenosos, a variação morfológica no aparelho de entrega de veneno, como presas e ferrões, pode ser otimizada para despachar presas. Aqui, determinamos como uma única dimensão dos sistemas de injeção de veneno evolui em resposta à variação no tamanho, condições climáticas e ecologia alimentar de cobras viperídeos. Medimos o comprimento das presas em mais de 1.900 espécimes de museus representando 199 espécies de víboras (55% das espécies reconhecidas). Encontramos tanto o sinal filogenético quanto a variação dentro do clado no comprimento relativo das presas entre as víboras, sugerindo tendências taxonômicas gerais e potencial divergência adaptativa no comprimento das presas. Recuperamos alometria evolutiva positiva e pouca alometria estática no comprimento das presas. Presas proporcionalmente mais longas evoluíram em espécies maiores, o que pode facilitar a injeção de veneno em presas mais volumosas. Por fim, aproveitamos os dados climáticos e de dieta para avaliar os correlatos globais do comprimento das presas. Descobrimos que os modelos de evolução do comprimento das presas são melhorados através da inclusão de temperatura e dieta, particularmente na medida em que as dietas são dietas pesadas para mamíferos. Essas descobertas demonstram como a variação adaptativa pode surgir entre os componentes de sistemas complexos de captura de presas.

A new species of Leptopelis (Anura, Arthroleptidae) from the south-eastern slope of the Ethiopian Highlands, with notes on the Leptopelis gramineus species complex and the revalidation of a previously synonymised species.

A new ground-dwelling species of treefrog in the genus Leptopelis is described from the Harenna Forest in south-eastern Ethiopia. The description is based on morphology and acoustics and is supported by molecular data. The new species has a small body size, and the digital discs on fingers and toes are significantly more conspicuous than in other semi-fossorial members of the L. gramineus complex. It occupies forest habitats at lower altitudes and is separated ecologically and geographically from high-altitude species of the complex. One of them, a parapatric cryptic species from Bale and Arsi Mountains, is resurrected from synonymy of L. gramineus and given a new name, L. montanus. Genetic barcoding of specimens from both populations showed that they belong to two distinct lineages that had been revealed by recent phylogenetic research. To confirm the geographic separation of the studied populations, the collection area of L. gramineus types was verified through analysis of the diary and the final report of the 2nd expedition of V. Bottego, and through matching of the route described in it with modern maps. The type locality of L. gramineus sensu stricto is restricted to Gamo Gofa, Ethiopia. Following the results of recent phylogenetic studies, the range of L. gramineus is limited to west of the Great Rift Valley. An identification key to the named Ethiopian species of the genus is provided.

PLA2 Inhibitor Varespladib as an Alternative to the Antivenom Treatment for Bites from Nikolsky's Viper Vipera berus nikolskii.

Although envenoming by a small East European species of viper is rarely severe, and only exceptionally fatal, lack of specific antivenom stocks in a few areas within this region and possible severe side effects of antivenom application leave most bites to be treated only with antihistamines and supportive therapy. Varespladib is an effective inhibitor of snake phospholipase, and, as such, it could be considered as first-line therapy. The Nikolsky's viper venom contains an extremely high concentration of phospholipase A2 (PLA2), responsible for the toxic effects of the venom, as well as minor amounts of other toxins. If Varespladib can successfully inhibit PLA2 activity, the Nikolsky's viper could be one of the first venomous snakes having an antitoxin-specific treatment regimen. To assess that, Varespladib was administered alone subcutaneously to adult male CD-1 mice (8 mg/kg) and compared to mice exposed to Vipera berus nikolskii crude venom (8 mg/kg = 10 LD50) or a combination of Varespladib and the same amount of the venom. Experimental animals were monitored for the presence of envenoming symptoms and mortality for 48 h after injection. Eighty percent of mice receiving both Varespladib and venom survived, while 100% of the control group receiving venom alone died within 4 h. Experimental results are consistent with Varespladib acting as an effective antitoxin in the mouse model against Nikolsky's viper venom. Further studies are needed under experimental conditions that more closely resemble natural envenoming (i.e., delayed administration).

Variation in body size and sexual size dimorphism in the most widely ranging lizard: testing the effects of reproductive mode and climate.

Reproductive mode, ancestry, and climate are hypothesized to determine body size variation in reptiles but their effects have rarely been estimated simultaneously, especially at the intraspecific level. The common lizard (Zootoca vivipara) occupies almost the entire Northern Eurasia and includes viviparous and oviparous lineages, thus representing an excellent model for such studies. Using body length data for >10,000 individuals from 72 geographically distinct populations over the species' range, we analyzed how sex-specific adult body size and sexual size dimorphism (SSD) is associated with reproductive mode, lineage identity, and several climatic variables. Variation in male size was low and poorly explained by our predictors. In contrast, female size and SSD varied considerably, demonstrating significant effects of reproductive mode and particularly seasonality. Populations of the western oviparous lineage (northern Spain, south-western France) exhibited a smaller female size and less female-biased SSD than those of the western viviparous (France to Eastern Europe) and the eastern viviparous (Eastern Europe to Far East) lineages; this pattern persisted even after controlling for climatic effects. The phenotypic response to seasonality was complex: across the lineages, as well as within the eastern viviparous lineage, female size and SSD increase with increasing seasonality, whereas the western viviparous lineage followed the opposing trends. Altogether, viviparous populations seem to follow a saw-tooth geographic cline, which might reflect the nonmonotonic relationship of body size at maturity in females with the length of activity season. This relationship is predicted to arise in perennial ectotherms as a response to environmental constraints caused by seasonality of growth and reproduction. The SSD allometry followed the converse of Rensch's rule, a rare pattern for amniotes. Our results provide the first evidence of opposing body size-climate relationships in intraspecific units.

Fifteen shades of green: The evolution of Bufotes toads revisited.

The radiation of Palearctic green toads (Bufotes) holds great potential to evaluate the role of hybridization in phylogeography at multiple stages along the speciation continuum. With fifteen species representing three ploidy levels, this model system is particularly attractive to examine the causes and consequences of allopolyploidization, a prevalent yet enigmatic pathway towards hybrid speciation. Despite substantial efforts, the evolutionary history of this species complex remains largely blurred by the lack of consistency among the corresponding literature. To get a fresh, comprehensive view on Bufotes phylogeography, here we combined genome-wide multilocus analyses (RAD-seq) with an extensive compilation of mitochondrial, genome size, niche modelling, distribution and phenotypic (bioacoustics, morphometrics, toxin composition) datasets, representing hundreds of populations throughout Eurasia. We provide a fully resolved nuclear phylogeny for Bufotes and highlight exceptional cyto-nuclear discordances characteristic of complete mtDNA replacement (in 20% of species), mitochondrial surfing during post-glacial expansions, and the formation of homoploid hybrid populations. Moreover, we traced the origin of several allopolyploids down to species level, showing that all were exclusively fathered by the West Himalayan B. latastii but mothered by several diploid forms inhabiting Central Asian lowlands, an asymmetry consistent with hypotheses on mate choice and Dobzhansky-Muller incompatibilities. Their intermediate call phenotypes potentially allowed for rapid reproductive isolation, while toxin compositions converged towards the ecologically-closest parent. Across the radiation, we pinpoint a stepwise progression of reproductive isolation through time, with a threshold below which hybridizability is irrespective of divergence (<6My), above which species barely admix and eventually evolve different mating calls (6-10My), or can successfully cross-breed through allopolyploidization (>15My). Finally, we clarified the taxonomy of Bufotes (including genetic analyses of type series) and formally described two new species, B. cypriensis sp. nov. (endemic to Cyprus) and B. perrini sp. nov. (endemic to Central Asia). Embracing the genomic age, our framework marks the advent of a new exciting era for evolutionary research in these iconic amphibians.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis).

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Hybrid origin of European Vipers (Vipera magnifica and Vipera orlovi) from the Caucasus determined using genomic scale DNA markers.

BACKGROUND: Studying patterns of introgression can illuminate the role of hybridization in speciation, and help guide decisions relevant to the conservation of rare taxa. Vipera magnifica and Vipera orlovi are small vipers that have high conservation status due to their rarity and restricted distributions in an area of the Caucasus region where two other related species are present - V. kaznakovi and V. renardi. Despite numerous observations of hybridization between different species of small vipers, and the potential of a hybrid origin for V. magnifica and V. orlovi based on their distribution with respect to V. kaznakovi and V. renardi, hypotheses of a hybrid origin have not been formally tested. Here we generate genomic-scale data by performing next generation sequencing of double digest restriction-site associated DNA libraries, and use these multilocus data to test whether these two species are of hybrid origin. RESULTS: We generated over nine hundred loci for 38 specimens of six taxa, and analysed the dataset using Bayesian clustering and multivariate methods, as well as Patterson D-statistics, which can distinguish between incomplete lineage sorting and introgression as explanations for shared polymorphism. The results demonstrate a pattern of historical admixture in the two purported hybrids that is consistent with past gene flow from V. renardi into V. kaznakovi. The average admixture proportion in individuals was low (6.39 %) in the case of V. magnifica, but was higher in V. orlovi (19.02 %). We also show that the specific individual samples used in D-statistic tests can have a significant impact on inferences regarding the magnitude of introgression, suggesting the importance of including multiple individuals in these analyses. CONCLUSIONS: Our results support the conclusion that both V. orlovi and V. magnifica had formed through a hybridization event between V. kaznakovi and V. renardi. Given a low proportion of admixture and absence of clear ecological and morphological differences V. magnifica should be treated as a marginal population of V. kaznakovi. Further studies that include analyses of ecological segregation of V. orlovi from parental taxa and search for evolutionary consequences of hybridisation would clarify if V. orlovi is a distinct hybrid species. Until this we recommend preserving the current taxonomy and protection status of V. orlovi.

Mitochondrial phylogeny shows multiple independent ecological transitions and northern dispersion despite of Pleistocene glaciations in meadow and steppe vipers (Vipera ursinii and Vipera renardi).

The phylogeny and historical demography of small Eurasian vipers of the Vipera ursinii and V. renardi complexes were studied using mitochondrial DNA sequences analysed with Bayesian inference, Maximum Likelihood and Maximum Parsimony approaches, and mismatch distributions. Diversification in the group resulted from an initial dispersion in the later Pliocene - Pleistocene in two directions: north-westwards via the Balkans (V. ursinii complex) and north-eastwards from Asia Minor via the Caucasus (V. renardi complex). An independent, comparatively recent transition occurred from montane habitats to lowland grasslands in different mitochondrial lineages during the Late Pleistocene, when representatives of the both complexes had reached lowland steppes to the north. Effective population size showed clear signs of rapid growth in eastern V. renardi, triggered by colonization of vast lowland steppes, but in western V. ursinii complex grew during the Last Glaciation and experienced stabilization in Holocene. Expansion and population growth in lowland lineages of V. renardi was not strongly affected by Pleistocene climatic oscillations, when cold, dry conditions could have favoured species living in open grasslands. The high diversity of closely related haplotypes in the Caucasus and Tien-Shan could have resulted from repetitive expansion-constriction-isolation events in montane regions during Pleistocene climate fluctuations. The mitochondrial phylogeny pattern conflicts with the current taxonomy.

Variation of Reproductive Traits and Female Body Size in the Most Widely-Ranging Terrestrial Reptile: Testing the Effects of Reproductive Mode, Lineage, and Climate.

The European common lizard, Zootoca vivipara, is the most widespread terrestrial reptile in the world. It occupies almost the entire Northern Eurasia and includes four viviparous and two oviparous lineages. We analysed how female snout-vent length (SVL), clutch size (CS), hatchling mass (HM), and relative clutch mass (RCM) is associated with the reproductive mode and climate throughout the species range and across the evolutionary lineages within Z. vivipara. The studied variables were scored for 1,280 females and over 3,000 hatchlings from 44 geographically distinct study samples. Across the species range, SVL of reproductive females tends to decrease in less continental climates, whereas CS corrected for female SVL and RCM tend to decrease in climates with cool summer. Both relationships are likely to indicate direct phenotypic responses to climate. For viviparous lineages, the pattern of co-variation between female SVL, CS and HM among populations is similar to that between individual females within populations. Consistent with the hypothesis that female reproductive output is constrained by her body volume, the oviparous clade with shortest retention of eggs in utero showed highest HM, the oviparous clade with longer egg retention showed lower HM, and clades with the longest egg retention (viviparous forms) had lowest HM. Viviparous populations exhibited distinctly lower HM than the other European lacertids of similar female SVL, many of them also displaying unusually high RCM. This pattern is consistent with Winkler and Wallin's model predicting a negative evolutionary link between the total reproductive investment and allocation to individual offspring.